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ABSTRACT 

The role of star-formation driven outflows in the obscuration of the central source 
in the Active Galactic Nuclei (AGN) is discussed. The outflow from a sub-parsec scale 
accretion disc is numerically modelled for parameters appropriate to the Galactic 
Centre. The resulting obscuration pattern is very patchy, with some lines of sight 
becoming optically thick to Thomson scattering. A flxed observer would see column 
depth changing by factors of many over time scales of order months to hundreds of 
years, depending on the physical size of the outflow region. Such winds may be relevant 
for obscuration of some AGN and especially "changing look AGN" . However, averaged 
over the sky as seen from the central source, these winds are always optically thin 
unless wind outflow rates are super-Eddington. A simple scaling argument shows that 
this is true not only for stellar-driven winds but for any AGN winds. We therefore 
conclude that AGN winds are unable to account for the vast majority of optically 
thick obscured AGN (a signiflcant fraction of all AGN). We suggest that the most 
likely source of optically thick obscuration in AGN is a warped parsec scale accretion 
disc. 

Key words: Galaxy: centre - accretion: accretion discs - galaxies: active - methods: 
numerical 



1 INTRODUCTION 

According to the simplest version of the Unified model of 
Active Galactic Nuclei (AGN) , the central source has funda- 
mentally the same properties in all classes of AGN (e.g., An- 
tonucci & Miller, 1985; Antonucci, 1993; Maiolino & Rieke, 
1995; Risaliti et al., 1999; Sazonov & Revnivtsev, 2004). It is 
the orientation of our line of sight with respect to the dusty 
obscuring environment of an AGN that makes it appear as a 
type I (relatively un-obscured) or type II object (obscured). 
Despite of all the convincing observational evidence for the 
existence of the absorber, most frequently thought to be ar- 
ranged in a shape of a torus, no convincing theoretical model 
has ever been produced to explain its properties. Recently, 
Nenkova et al. (2002) and Risaliti et al. (2002) demonstrated 
that to comply with the small observed torus sizes, the ab- 
sorber must not be uniform in density but should be rather 
clumpy, i.e. consisting of individual smaller clouds. Theoret- 
ical models of clumpy tori (e.g., Krolik & Begelman, 1986, 
1988; Vollmer et al., 2004), however, have all the same fun- 
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damental difficulty. As the observations require the torus to 
be geometrically thick, the velocity dispersion of the cold 
dusty clouds must be of the order of their orbital velocity 
in the inner parsec of a galactic centre. Optically thick tori 
would then have clouds colliding at super-sonic speeds many 
times during one dynamical time (Krolik & Begelman, 1986, 
1988). Therefore, these models have to postulate (somewhat 
arbitrarily, in our point of view) , that the clouds can survive 
these collisions, behaving in essence as rubber balls. 

Recently, a number of authors argued that the obscu- 
ration must instead come from clumpy winds (e.g., Konigl 
& Kartje, 1994; Kartje et al., 1999; Elvis, 2000; Elvis et al., 
2004; Elitzur, 2005) emanating from the surface of an ac- 
cretion disk, be these winds driven by line pressure (Proga, 
2003), continuum radiation pressure, heating, or by hydro- 
magnetic forces. In this paper we wish to critically assess this 
suggestion. In particular, we consider AGN winds driven by 
star formation activity in an accretion disc. Star or planet 
formation in AGN discs has been long predicted by theo- 
rists (Paczynski, 1978; Kolykhalov & Sunyaev, 1980; Lin & 
Pringle, 1987; Collin & Zahn, 1999; Gammie, 2001; Good- 
man, 2003), and recent observations of the young massive 
stars in the centre of our Galaxy (e.g., Genzel et al., 2003; 



2 S. Nayakshin and J. Cuadra 



Ghez et al., 2003) lend a very strong support to these the- 
ories (Levin & Beloborodov, 2003; Nayakshin & Sunyaev, 
2005; Paumard et al., 2006). Using the SPH treatment of 
stellar winds developed by Cuadra et al. (2005, 2006), we 
calculate the obscuring properties of such winds for a spe- 
cific example. We find that these winds have very irregular, 
patchy obscuration patterns. A fixed observer should then 
see the column depth through the wind varying on time 
scales from months or years to tens and hundreds of years, 
depending on the physical size of the wind launching region 
and the SMBH mass. 

However, observations (see §3 below) suggest that the 
structure responsible for AGN obscuration should be opti- 
cally thick over a large (e.g., a third) fraction of the whole 
sky. Our models cannot achieve such a high average line of 
sight column depth unless the mass loss rate in the wind is 
super-Eddington. In fact, we show via simple analytical esti- 
mates that this conclusion is valid whatever the nature of the 
wind driving mechanism is. We feel it would be very hard to 
produce highly super-Eddington winds in a plausible quasi 
steady-state AGN model for a typical (very sub-Eddington 
luminosity) AGN. Thus we suggest that winds, while very 
important for the AGN phenomenon, cannot replace the 
"torus" in its role in the unification schemes of AGN. Ei- 
ther the torus does exist despite the difficulties faced by the 
models, or another optically thick structure, e.g., a warped 
accretion disc plays its role in AGN. 



2 STELLAR WINDS FROM AGN DISCS 

AGN discs are expected to be massive and thus self- 
gravitating at large enough distances (i?^ 0.01 — 0.1 pc) 
from the super-massive black holes (SMBH). If the disc 
cooling time is short, gravitational collapse and formation 
of stars or even planets is predicted (e.g., Paczyhski, 1978; 
Kolykhalov & Sunyaev, 1980; Shlosman & Begelman, 1989; 
Gammie, 2001). Young massive stars in the central parsec 
of our own Galaxy were most likely created in this way 
(Paumard et al., 2006), with an apparent significant over- 
abundance of high mass stars (Nayakshin & Sunyaev, 2005) 
over their fraction in "normal galactic" star formation. As 
stars are born from the accretion disc, they launch powerful 
radiation fields and stellar winds, some of which will break 
through the disc. Low mass proto- stars may prove equally 
effective in launching winds in these circumstances as the 
rates at which gas is captured from the disk into the Hill 
(or capture) zones of the protostars are super-Eddington 
(Nayakshin, 2006). In addition, when the gas in the disk 
is depleted to about half its initial surface density, stellar 
velocity dispersion starts to grow, and soon the stellar disc 
becomes geometrically thicker than the gas disc (Nayakshin, 
2006). Stellar winds then escape from above the disc directly. 

We attempt to simulate this complex situation in a sim- 
pliffed setting, concentrating only on the wind part of it. For 
this reason we do not include the gas in the negligibly thin 
and ffat accretion disc from which the stars are born. To 
capture a degree of the expected diversity of the stellar pop- 
ulations and conditions in this problem, the stars are divided 
into two groups. The ffrst produces winds with velocities 
I'w = 300 km/sec, whereas the second has v^^ = 700 km/sec. 
Both types of stars have mass loss rates of 2.5 x 10~^ M© 



year~^, and are situated in a ffat circularly rotating Kep- 
lerian disk of geometrical thickness H(R) — O.IR (the unit 
of length used here is 1 arcsecond at the distance of 8 kpc, 
which is about 1.2 x 10^^ cm or 0.04 pc). In total, we have 200 
mass shredding stars, thus amounting to the mass loss rate 
of 0.1 Mq year~^, which is a factor of few super-Eddington 
for a SMBH mass of Mbh = 3.5 x 10^ M©. The disc inner 
and outer radii are Rin = 1.5 and i^out = 8, respectively. 
The stellar surface density follows the law S*(i?) oc 
and the stars are in the Keplerian circular rotation around 
the SMBH. The initial phase (i.e. the 0-coordinate) of the 
stars in the disk is generated randomly, as is the vertical co- 
ordinate within the —H to -\-H limits. Initially, the compu- 
tational domain is filled with hot tenuous gas with velocity 
greater than the escape velocity. This gas quickly outfiows 
from the region. At the end of the simulation stellar winds 
exceed the initial gas mass by a factor of few tens. 

To perform simulations, we use the Gadget- 2 code 
(Springel, 2005) modified to model stellar winds as new SPH 
particles ejected from the stellar wind sources. These meth- 
ods were described and thoroughly tested in Cuadra et al. 
(2005, 2006). 

A snapshot of the simulation is shown in Figure 1 at 
time t ~ 2200 years after the beginning of the simulation, 
which corresponds to ~ 36 dynamical times at R = 1. While 
there is no true steady state in this finite number of moving 
stars system, the snapshot is fairly typical of the morphology 
of the stellar wind. The face-on view (left panel) shows that 
some of the shocked wind managed to cool down and formed 
a small-scale disc. The inner part of the simulation domain is 
the most likely place for the disc to form since stellar wind 
density is highest there, leading to shocks, thermalization 
and rapid cooling; the frequent collisions of gaseous clumps 
"waste" momentum of the winds, and finally, the escape 
velocity from the region is around 600 km/sec, i.e. higher 
than wind velocities of the slower winds. With time the disc 
grows radially to both smaller and larger radii. However, the 
formation of the disc (cf. Cuadra et al., 2006) is physically 
significant only in the case when the gaseous accretion disc 
from which the stars were born is entirely absent, as is the 
case presently in the Galactic Centre. In the opposite case 
the cooled disc and the "spiral arms" seen at intermediate 
radii would simply blend in with the much more massive 
underlying gas accretion disc. 

Contrary to the frequent collisions and high escape ve- 
locity in the inner region of the computational domain, the 
conditions in the outer regions allow direct escape of both 
the faster diffuse and the slower cooler clumpy winds. Due 
to the final extent of the stellar disc and the projection ef- 
fect, the wind morphology reminds and "X" -shape (Figure 
1, right panel). The edge-on view (right panel) of the stellar 
disc shows that the plane of the gas disc at ^ 2 is very 
slightly tilted with respect to the orbital plane of the stellar 
disc. This is an artefact of the initial conditions. 

The left panel of Figure 2 shows the obscuring column 
depth of the winds as seen from the SMBH. Since the orien- 
tation of the tilted inner disc is infiuenced by the initial con- 
ditions, we eliminate its obscuration, only plotting gas with 
radial distances ^ 1.6 from the SMBH, which excludes all 
of the disc. Besides, had we included in the simulations the 
more massive gas accretion disk whose midplane would co- 
incide with that of the stars, the inner disc seen in Figure 
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Figure 1. Face-on (left) and side (right) views of the simulation domain for the flat stellar disc conflguration. Red asterisks show the 
location of the wind-producing stars used in the simulation (the red "dot" in the centre of the left panel is not a star, but a low column 
depth region centred on the SMBH). 



1 would have completely blended in with the accretion disc 
and would not present much obscuration at all. 

Note the very irregular patchy structure of the (yel- 
low) optically thicker regions. The contrast between those 
and neighbouring less dense patches of sky is frequently 
a factor of 10 or more. Since the winds are rotating at 
a fraction of the local angular frequency, Q, this implies 
that the column depth sampled by the observer will vary 
on time scales as short as ^ 10~^ to 10~^ of l/H, i.e. 
tvar 0.04 - 0.4 year {R/O.lpcf^'^ M~^^^ . The dotted pat- 
tern at the cos^ = plane are the dense regions of stellar 
winds immediately next to the stars. 

The right panel of Figure 2 shows the isotropic mass 
loss rate along the line of sight, defined as 



M : 



7^^ 



(1) 



where the integral is taken along the line of sight defined 
by a given and cj). This quantity is useful for comparison 
with observations of AGN outflows. An observer has only 
one line of sight available at any given moment, and the 
assumption frequently made is that the outflow is roughly 
isotropic, thus M ~ AttR^pvr. We find that there are lines of 
sight that yield isotropic mass outflow rates of '-^ 1 — 10 M© 
year~^, which is as much as a hundred times larger than 
the correct sky- averaged value. Also note that not all of the 
optically thicker regions seen in the left panel of figure 2 
appear equally prominently in the outflow map in the right 
panel of the figure, as some of these structures have a small 
radial velocity or are even infalling to smaller radii. 



3 COMPARISON TO OBSERVATIONS 

Recent multi- wavelength observations (e.g., Jaffe et al., 
2004; Packham et al., 2005; Prieto et al., 2005) indicate that 
AGN absorbers are relatively small, parsec-scale structures, 
rather than the extended (100 to 1000 parsec) ones as was 
thought a decade earlier on the basis of observations with 
poorer angular resolution. Nenkova et al. (2002) and Elitzur 
(2005) demonstrate that this observational fact demands the 
"torus" to be clumpy in order to allow high and low tem- 
perature regions to be simultaneously present at the same 
location in the torus. Nenkova et al. (2002) (see also Elitzur, 
2005) show that the average number of clouds on the Une 
of sight to the nucleus, A^a, should be of the order of a few 
to ten. Let Mc be the mass of such a cloud, Mc = Sc7ri?c, 
and Ec and Rc be the surface density and cloud radius, re- 
spectively. In the model of these authors, the average column 
depth on the line of sight, N^^c is of the order of 1 g cm~^ or 
larger, which is also reasonable given a large fraction of AGN 
that are Thomson-thick (e.g., Sazonov & Revnivtsev, 2004; 
Guainazzi et al., 2005). The total number of the clouds. At, 
is connected to A^a through 



ATa^ 



NtTTRj 

21TRI 



(2) 



where Rt is the typical radius of the torus, and we assumed 
that it covers roughly 2n of the sky as seen from the SMBH. 
Now, assume that the torus is made up of clouds that form 
in the outflowing disk wind, as in the simulations presented 
in this paper. The total mass outflow rate of the clouds is 
then 
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Figure 2. Column depth through the wind (left) and isotropic mass outflow rate (right; in units of year"-*^) for the simulation 
shown in Fig. 1. Note the large variations of the obscuring column depth over small angular scales. The gaseous inner disc seen in the 
left panel of Figure 1 has been excluded from these Figures. 



Mwind - 27vR^,N,^c^K - 15^ iVa5]c ^ Mg'/' ,(3) 

year yl pcy 

where Ms is the SMBH mass in units of 10^ M©, and we 
assumed that clouds radial velocities are of the order of 
the local Keplerian speed. Now, the Eddington accretion 
rate is AttG MBurrip / ecar ~ 2M0year~^M8. Thus the re- 
quired wind mass loss rate (equation 3) is an order of mag- 
nitude higher, typically, than the Eddington accretion rate 
for a Ms ^ 1 object and a parsec-scale torus. The estimate 
of Mwind could be reduced somewhat by postulating even 
smaller torus sizes, but it is hard to make Rt much smaller 
than 0.1 pc as then the dust would be sublimated by the 
AGN radiation field (e.g., eq. 3.2 in Emmering et al., 1992). 
However, we rather think that estimate 3 is too optimistic, 
e.g. that an even higher mass outflow rates are needed, as 
the observations and modelling of optically thick sources re- 
quire a minimum absorbing column depth, and hence in a 
good fraction of sources N^T^c may actually be much higher 
than 1 g cm~^. 

Considering a specific case of the local obscured AGN 
studied by (Guainazzi et al., 2005), we note that the bolo- 
metric luminosities of these objects in the infrared. X-ray 
and optical bands are in the range L ^ 10^^ — few x lO'^^ 
erg/sec, which implies SMBH accretion rates of "only" ^ 
0.01 M0 year~^ for the standard radiative efficiency. Hence 
if the obscuration of the optically thick objects in that sam- 
ple were provided by the winds, we would conclude that 
the SMBH accretion process must be very wasteful, with 
~ 100 — 10, 000 times more mass flowing out of the inner 
parsec than accreting on the SMBH. (Note that these mod- 
erately bright AGN are not likely to be in the non-radiative 



accretion flow regime when vigorous outflows are in fact ex- 
pected, e.g., Blandford & Begelman, 1999) It would also re- 
quire a very high mass influx into the inner parsec to sustain 
such winds. Given the difficulty of delivering enough fuel to 
the SMBHs even in the earlier gas-rich epochs (Thompson 
et al., 2005), it is hard to see how such high mass influxes 
could be maintained. 



4 DISCUSSION AND CONCLUSIONS 

Accreting black holes can drive strong outflows via X- 
ray heating, line or continuous radiation pressure, and hy- 
dromagnetic forces (e.g., Begelman et al., 1983; Konigl & 
Kartje, 1994; Kartje et al., 1999; Proga, 2003). Star forma- 
tion on the outskirts of cool massive accretion discs will also 
result in winds driven by outflows from the young massive 
stars (Cuadra et al., 2005). Here we studied the obscura- 
tion properties of the AGN outflows of the last type. We 
found that these outflows are quite clumpy, with the lines 
of sight passing through the clumps becoming moderately 
optically thick to Thomson scattering for high enough wind 
mass loss rates. Such outflows are bound to play a role in 
the obscuration of AGN. 

We however feel that this role cannot be dominant 
as the sky averaged column depth is significantly smaller 
than the one required by the current "torus" models 
(e.g., Nenkova et al., 2002; Elitzur, 2005) and observations 
(Guainazzi et al., 2005) for realistic mass outflow rates, e.g. 
comparable to the accretion rates onto the SMBH, which 
are likely to be only 0.01 — 0.1 of the Eddington accretion 
rates in the local AGN (and are much smaller if we also 
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Figure 3. Same as the left panel of Figure 2 but for slightly 
different initial conditions and with the cold tilted and warped 
inner disc included in the Figure. Notice that this very light (only 
~ 3M0) disc is a very efficient absorber. 



consider Low Luminosity AGN; see, for example, Fig. 1.9 in 
Ho, 2005). We thus believe that the search for the culprit of 
the optically thick obscuration of AGN should still go on. 

One rather obvious candidate is a warped accretion disc. 
The disc may be significantly warped by instabilities due to 
back reaction to mass outflow or due to the AGN radiation 
pressure, or it may be warped due to precession in a non- 
axisymmetric gravitational potential (e.g. Schandl & Meyer, 
1994; Pringle, 1996; Nayakshin, 2005). Finally, there is no 
reason for an initial disc configuration to be flat if the mass 
deposition is much shorter than the disc viscous time, which 
was probably the case in the Galactic Centre (Nayakshin & 
Cuadra, 2005; Paumard et al., 2006). A simulation ran by us 
for the present paper and later rejected (because the initial 
mass of the gas turned out too high so that it influenced the 
outcome of the simulation) is curiously useful in this respect. 
Figure 3 shows the obscuration pattern of this simulation in 
the same way as in the left panel in Figure 2, but including 
the tilted and slightly warped gas disc. Except for the initial 
conditions, this simulation is exactly identical to the one 
shown in Figures 1 & 2. As Figure 3 shows, the obscuration 
provided by the winds literally pales in comparison with 
that produced by the warped disc, whose mass is only ^ 
3 M0 . An accretion disc is a long lived structure that can 
accumulate a lot of gas (still a small fraction of the SMBH 
mass, however) and can be Thomson-thick out to parsec 
distances. Obscuration by a warped accretion disc requires 
no continuous energy input, and would work for low or high 
luminosity AGN as long as the discs are sufficiently massive 
and strongly warped. 
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